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Abstract 

Signatures of glassy dynamics have been identified experimentally for a rich variety of materials 
in which molecular networks provide rigidity Here we present a theoretical framework to study the 
glassy behavior of both passive and active network materials. We construct a general microscopic 
network model that incorporates nonlinear elasticity of individual filaments and steric constraints 
due to crowding. Based on constructive analogies between structural glass forming liquids and 
random field Ising magnets implemented using a heterogeneous self-consistent phonon method, 
our scheme provides a microscopic approach to determine the mismatch surface tension and the 
configurational entropy which compete in determining the barrier for structural rearrangements 
within the random first order transition theory of escape from a local energy minimum. The 
influence of crosslinking on the fragility of inorganic network glass formers is recapitulated by the 
model. For active network materials, the mapping, which correlates the glassy characteristics to the 
network architecture and properties of nonequilibrium motor processes, is shown to capture several 
key experimental observations on the cytoskeleton of living cells: Highly connected tense networks 
behave as strong glass formers; intense motor action promotes reconfiguration. The fact that our 
model assuming a negative motor susceptibility predicts the latter suggests that on average the 
motorized processes in living cells do resist the imposed mechanical load. Our calculations also 
identify a spinodal point where simultaneously the mismatch penalty vanishes and the mechanical 
stability of amorphous packing disappears. 
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I. INTRODUCTION 



Materials in which molecular networks provide rigidity are ubiquitous, ranging from pas- 
sive systems such as rubber and silica glasses to active systems, such as the crosslinked 
polymer networks driven by energy-consuming motor proteins that constitute the cytoskele- 
ton of eukaryotic cells. Both active and passive network materials exhibit features of glassy 
dynamics in parts of their phase diagrams. Eukaryotic cells possess a cytoskeleton that 
stiffens under tension while having an intracellular space crowded with macromolecules and 
organelles that can resist compression. Recent experiments have established striking simi- 
larities between the behavior of a living cytoskeleton and that of inert nonequilibrium soft 
glasses^. Studies of the mechanical and dynamical behavior of osmotically compressed cells 
further suggest that cells under compressive stress behave as strong glass formers^. At 
the same time, F-actin disruption and ATP-dependent nonequilibrium processes have been 
shown to strongly modulate the glass transition behavior. 

Here we develop a framework for studying the glassy behavior of network materials and 
report calculations for a model of the cytoskeleton to illustrate the concepts. Our scheme is 
based on the explicit analogy between a structural glass forming liquid and a short range dis- 
ordered Ising magnet used earlier to study the Lennard- Jones glass^. This mapping correlates 
the glassy characteristics with the network architecture (density and connectivity) and the 
properties of motor-driven active processes. The approach provides an explicit microscopic 
route to calculate from the force laws the configurational entropy and the mismatch energy 
or surface tension, which are competing factors in determining the barrier for structural 
rearrangements within the random first order transition (RFOT) theory^. Site-dependent 
Debye- Waller factors, which characterize the extent of localized motion, are determined by 
a heterogeneous self-consistent phonon (SCP) method. Our calculations identify a spinodal 
point at a dynamic transition packing fraction where the mismatch surface tension and 
mechanical stability simultaneously vanish^"—. 

Combining the magnetic analogy and self-consistent phonon method allows us to build 
a correspondence between the theory and several experimental observations: The statistics 
of the random fields and random interactions calculated from the present network model 
indicates a possible ideal glass transition underlying the observed glassy behavior of highly 
connected tense networks made of stiff filaments. When there is a higher degree of connec- 
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tivity a crossover to strong liquid behavior takes place, both for passive and active materials. 
Load-resisting motor processes are shown to facilitate structural rearrangements by expand- 
ing the accessible configurational state space and lowering the reconfiguration barrier. 

II. A GENERAL MICROSCOPIC NETWORK MODEL 

The microscopic model we use can be thought of as an amorphous "cat's cradle" with 
excluded volume^^^— . It consists of a crosslinked network of nonlinear elastic filaments that 
stretch elastically with effective stiffness /37 beyond the relaxed length Lg but that will 
buckle and bear no load if shortened too much. We model the crosslinks as hard spheres 
of diameter dns that provide the excluded volume. In the cytoskeletal instantiation this 
excluded volume represents the steric constraints from both the filaments and the binding 
proteins. Instead of following the filament degrees of freedom, we keep track of the motion of 
the crosslinks or the nodes of the network. Bonded node pairs interact with the cat's cradle 
type potential given by (5U{r) = B(r — Le)(l/2)/37(r — Lg)^ + AQ{dHs — '^)) where B(-) is 
the Heaviside step function and A ^ oo indicates hard-core repulsion. Elastic energy arises 
only when the contour length r exceeds the relaxed length Lg. The statistical architecture of 
the network is specified by the node density po and network connectivity (with Pc denoting 
the fraction of bonded neighbors). We assume the network architecture is quenched once 
initially assigned. 

To deal with active network materials we have shown earlier^^ that when agitated by 
spatially anti-correlated motor-driven events the cytoskeleton behaves as if it were at an 
effective equilibrium with a nontrivial effective temperature and modified interactions. The 
effective motor-driven pair interaction scaled by the effective temperature Teff (= is 
given by /3efft4fr(^) = [1 + (^^ ~ V^)^] PU{r) + Alnr, which depends on motor activity A 
and susceptibility s. A schematic illustration of the pair interaction in passive and motor- 
driven active networks is shown in Fig. [H In contrast to the passive network (blue curve), 
motor action induces an effective attraction (A Inr term) even in the buckling regime {dns < 
r < shaded area). Motor susceptibility modulates the long-range elastic interactions: 
Susceptible motors with 5 > (purple curve) enhance long-range attraction, whereas intense 
action of load-resisting motors with 5 < (magenta curve) may yield long-range repulsion. 

Using the underlying parameters of the model and starting from an arbitrary initial con- 
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FIG. 1. Schematic of model interactions. Blue line: pU{r) for passive networks. Elastic stretching 
initiates at r = Lg. Purple (s > 0) and Magenta (s < 0) lines: /3eff^^efr(^) foi" active networks. 
Motor action induces effective attraction even in the buckling regime (shaded area). Susceptible 
motors (s > 0) enhance the long range attraction whereas load-resisting motors {s < 0) may lead 
to long-range repulsion at high motor activity. 

figuration, one can perform molecular dynamics (MD) simulations to obtain representative 
fiducial structures for a given set of parameters (po, ^e, /^t; A, 5). These structures are 
vibrationally equilibrated but quenched as to density and connectivity. We fix the number 
density po = 1-2 of a total of = 256 particles, setting the mean particle spacing tq as the 
length unit. Periodic boundary conditions are applied. 

We start with a disordered configuration from a simulated Lennard- Jones binary mixture 
used for earlier studies of glassy dynamics^ but assume a single radius in the starting config- 
uration for current purposes. We then evolve the configuration via MD steps until the mean 
square displacement (MSD) of the nodes with respect to their initial positions saturates to 
a plateau (see an example in Fig. [2^). During equilibration a series of jerks in the MSD 
arise in the absence of thermal noise (zero temperature assumed). These may be thought 
of as avalanches, corresponding to large-scale cooperative rearrangements. During the MD 
evolution toward the steady state structure, as seen in Fig. [2)3, the repulsive interaction 
energy (Er) decreases with time whereas the attractive interaction energy (Ea) due to bond 
stretching increases as structural rearrangements occur; the total potential energy (Efot) 
decreases as time advances until an apparent steady state plateau is reached. The radial 
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distribution functions g{r) for the initial and final configurations are given in Fig.[2]3. These 
show how the original binary sphere system (red line) with two peaks in the innermost shell 
evolves into a monodisperse network structure (blue line) with a single peak in the first shell. 
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FIG. 2. Monitoring the generation of fiducial structures, (a) Mean square displacement (MSD) 
versus number of MD steps, (b) Evolution of the potential energy. Ea and denote for attractive 
and repulsive interactions, respectively, and Etot for their sum. (c) Radial distribution function 
g(r). Red: initial binary system; blue: final monodisperse structure. Model parameters are Lg = 
1.2,/37 = 2.5,77 = 1. 



III. MAGNETIC ANALOGIES OF THE RFOT THEORY 



A structural glass is statistically homogeneous but is specifically nonuniform at the mi- 
croscopic level with a density p{x) that is not translationally invariant. This is what gives 
rise to the glass's rigidity. Liquid state theory provides the free energy as a function of 
such a non-uniform density F[p{x)]. The one-to-one correspondence in the long time limit 
between force field and density profile has been elegantly established a la the density func- 
tional theory—. While the complete equilibrium free energy assumes the whole phase space 
can be sampled, close to the glass transition there occurs a trapping in locally metastable 
states which manifest themselves as an extensively large number of local free energy minima 
described by a configurational entropy Sc = InA^^^^. Self-consistent phonon (SCP) theory 
directly addresses broken ergodicity by investigating the stability of a density wave frozen 
into aperiodic structures. 
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The random first order transition (RFOT) theory relates the activation energy of a rear- 
ranging unit (the "droplet") to the microscopic intermolecular forces, and has successfully 
predicted many confirmed quantitative results for the dynamics of glass forming liquids^. 
The essence of the RFOT theory of escape from a local minimum is an analogy to a random 
field Ising magnet in a biasing field. Stevenson et al^ have established an explicit mapping 
between a structural glass forming liquid and the corresponding short range disordered Ising 
magnet. This mapping yields a description of a structural glass in terms of discrete spin-like 
variables tied to the liquid structure equilibrated at one time. The average field in the Ising 
magnet corresponds to the configurational entropy density Sc for the model liquid, which pro- 
vides an entropic driving force for the transitions among distinct configurational states. The 
smallness of Sc indicates a paucity of accessible configurations, and thus signals a deep de- 
scent into the glassy regime. Fluctuations in configurational entropy, Ssc^ on the other hand, 
are related to the excess heat capacity ACp (according to Landau, 5sc = ^/ /S^C^kB/Ncorr 
where Ncorr is the volume within which the disorder is correlated) and thus measures the 
fragility of the liquid. The strength of the coupling Jij between neighboring Ising spins that 
point in opposite directions, is associated with the interaction energy between a low-overlap 
site and its high-overlap neighbor or vice versa. 

In a three dimensional frozen aperiodic state, the density profile is well approximated 
by a sum of Gaussians centered around the fiducial lattice sites, p(r) = ^iPi{r) = 
'^.{ai/Tr)^^'^ exp[—ai{r — r-Y]^ where {a^} represent the effective local spring constants 
that determine the mean square displacement from the fiducial sites {r-}. For large values 
of {ai} the particles are localized close to the fiducial locations. By using the independent 
oscillator approximation which yields site- wise decoupling of the particles, the free energy 
^Fgiass can be expressed as a sum of effective potentials between the interacting density 
clouds, /3Veff{\f/ — fj\]aj) = — In J drjpj{rj)e~^^^^^i^~^^\ Within the Gaussian density 
ansatz one finds the free energy for the glassy state 

/3F,iU{rIh {«!}) = E 1 1^ ^ + E /^^^ffd^t - n'l; 4)- (1) 

J 1 

i ij 

Near the low overlap state, the free energy follows from the equilibrium liquid equation 
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of state (given by Zeos{v)) with corrections due to bonded constraints/interaction, 




Note the double summation in the last term only involves bonded pairs. 

For any assignment of the discrete values of {aj} and {af}^ the free energy of the model 
liquid is equivalent to a pairwise interacting model with spins located at the fiducial lattice 
sites, {f/}, 

i i<j 

where the spin Si = 1 corresponds to a large overlap site and 5^ = a small overlap 
site. The average field is found from the bulk free energy diflFerence between the states, 
"^ihi = /3Fgiass — f3Fiiq = Nsc/kB-) with a heterogeneous local field resulting from spatial 
variations of a, 

= I In 4 + I + E /^^effdrY - r^^l; 4) - fiZEosi^') - 1]^ 

- E[/5Kr^' - /5Kff^](l Y - rl\;a]). (4) 

j 

Here the summation in the last term, i.e. the modification due to bonding, includes only 
bonded neighbors of the central node i. 

The quenched interactions between neighboring low-overlap and high-overlap sites defined 
through the effective potential give the surface energies of the droplets within the RFOT 
picture and are explicitly 

J^J = PVf{\r/ - i\;ata])+pvf{\f/ - fl\;aj,aj), (5) 

where the pair interaction is given by 

Py2^{\r/ - rj\;al a]) = J dudfjpj{u)p]{fj)e-'^^^^^^^ (6) 

The statistics of the mapped out random fields (Eq. Hj) and random interactions (Eq. [5j), 
which depend on the liquid structure, when combined with the renormalization group results 
for random magnets^^ allow us to predict whether a thermodynamic phase transition would 
occur for the liquid analog when its mean field configurational entropy vanishes. 
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IV. HETEROGENEOUS SELF-CONSISTENT PHONON (SCP) METHOD 



The site-dependent Debye- Waller factors {a^} can be determined using the heterogeneous 
SCP method applied to the frozen/quenched fiducial lattices generated by the molecular dy- 
namics simulations. The validity of the mapping relies on the self-consistency between the 
fiducial structure and the inter-particle potential, in other words, the one-to-one correspon- 
dence between the density profile and the force field, as well as a well separation of time 
scales between the vibrational motion about a fiducial configuration and the collective struc- 
tural rearrangements representing hopping between different configurational states. Starting 
from uniform low (of order 1) and high (of order 100) a values, we obtain two distinct sets 
of mechanically stable homogeneous solutions {af} and {o^J}, stabilized by bond stretching 
and hard-core repulsion, respectively. An example of the color-coded a configurations is 
shown in the left two panels of Fig. [31 The modest spatial variation in a values that is seen 
refiects the disorder in quenched lattice structures and network connectivity. 

The SCP method also allows us to study situations where the order parameters spatially 
vary across the sample. In particular, we can study the case where a smooth interface forms 
between liquid-like and solid-like bulk stable phases. Fig. [3] illustrates the formation of a 
broad interface (see configuration marked as "final") from a mechanical relaxation (using 
SCP procedure) across an initially sharp interface (see configuration marked as "initial") 
between distinct mobile (reddish) and localized (bluish) bulk states. The a values within 
the left and right boundary layers (indicated by anchored arrows) have been pinned during 
the relaxation to maintain the spatial gradient in a. 

V. MAIN FINDINGS 
A. Passive networks 

1. Implications of the mapping 

Highly connected networks exhibit strong liquid behavior and descend deeper into the 
glassy regime 

We calculate the local fields {hi} defined by Eq. [4] using the homogeneous solutions {aj} 
and {aj} obtained by site-dependent SCP calculation. According to the mapping, h is 
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FIG. 3. Illustration of interface formation. A quenched typical fiducial particle configuration with 
color-coded a values is shown. Color scale is given on the right hand side. Green frames indicate 
the creation of an initial sharp contact between two bulk stable phases. With the a values in the 
boundary layers being pinned (indicated by blue and red anchored arrows), the final a configuration 
develops a smooth interface, showing a gradual progression of phases. 




FIG. 4. Highly connected networks descend deeper into the glassy regime and exhibit strong-liquid 
behavior, (a) The configurational entropy density Sc vs connectivity Pc- (b) The fluctuations 
of configurational entropy density Ss^ vs connectivity Pc- Model parameters are Lg = 1.2, /37 = 
2.5,7/ = 1. 

translated into the configurational entropy density, i.e., h = Sc/ks- In RFOT theory the 
configurational entropy parametrizes a liquid's descent into the glassy regime. We show in 
Fig.[4](a) that increasing Pc for a fixed density yields a smaller 5c, suggesting that crosslinking 
networks to a higher degree of connectivity leads them deeper into the glassy regime. This 
finding is consistent with earlier results using a statistically averaged RFOT-SCP theory 
for network glasses^^. The result also agrees with the well-known experimental observation 
that the addition of nonbonding impurities into highly crosslinked network glasses lowers the 
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glass transition temperature or raises the transition densities. Physically, higher connectivity 
yields a much stronger localization even for the loosely tethered (liquid-like) state; the 
resulting change in the entropy cost to localize density waves dominates over the reduction 
of bonded interaction due to smaller fluctuations in motion, leading to a higher liquid-state 
free energy and hence a lower Sc for the same conflguration. This is because the glassy-state 
free energy itself is insensitive to Pc since the caging of the motion arises primarily from 
steric constraints rather than bond stretching. 

On the other hand, the fluctuations of the conflgurational entropy, Ss^^ can be directly 
related to the excess heat capacity ACp. The mapping thus establishes the connection 
between a liquid's fragility and its degree of bonding: As seen in Fig. [4](b), higher Pc yields 
smaller 5sl and thus smaller AC^, indicating that liquids with a higher degree of bonding 
are stronger. Again this coincides with experience for inorganic network glass formers, in 
addition, this trend is consistent with the experimental observation that F-actin disruption 
in the cytoskeleton substantially increases its fragility-. 

^ Quasi-universality with respect to network connectivity in deeply glassy regime 

We have examined the statistics of the mapped out random interactions according to 
Eq. [51 Fig. [5](a) displays the distributions, at various levels of crosslinking, of the 
calculated interactions as the free energy per neighbor, Ji = — • Ja. where Zi is the co- 
ordination number of particle i. The typical interaction J = -^r Ji is related to the 
mismatch free energy penalty in RFOT theory for a particle at a flat interface between 
regions of high and low overlap. The Xia-Wolynes estimate^ for surface tension gives 
Jrfot = bT In -^^/n^i) 0.58, where rf^ = (l/10)^o is the Lindemann length and 
77,55 = 3.2 is the typical number of bonds broken by the interface. As shown in Fig. [5](b), 
the calculated J for single-particle droplets (low-a sites against a high-a glassy background) 
is insensitive to Pc and has a similar value to the Xia-Wolynes estimate. The absence of 
mechanical relaxation across the sharp interface would explain the modest overestimation of 
the energetic penalty. Fig. [5](c) demonstrates the quasi-universal proflle of the broad inter- 
face with respect to Pc if mechanical relaxation between the pinned boundaries is allowed, 
such that a gradual progression of phases lowers the energetic cost of forming an interface. 
Locating analogous magnets on the Renormalization Group phase diagram 

Analysis of the statistics of the calculated random flelds and interactions also allows us to 
locate the analogous magnets on the random fleld Ising magnet phase diagram determined 

10 




FIG. 5. Quasi-universality of the surface energy and the interface width with respect to network 
connectivity, (a) Distribution of the calculated local interactions {Ji} of the analogous magnets. 
Left to right: Pc = 0.5, 0.8 and 1. (b) J vs Pc- Typical surface tension J (= (1/A^) J2i Ji) dose to 
the RFOT estimate Jrfot — 0.58. (c) Initial and final a profiles at various A broad interface 
forms between localized and mobile bulk states. Model parameters are Pg — 1-2, /37 = 2.5,77 = 1. 

by renormalization group analysis at zero average field (Fig. [6j). Both the field fluctuations, 
and the fluctuations of the interaction strength, (5Jz^/^, are normalized by the total 
interaction energy per site Jz. We use the Ji distributions determined for single-particle 
droplets. 

In Fig.Elthe colored symbols indicate estimates of where the analogous magnets for model 
networks with various nonlinearity parameters and connectivity would fall on the magnetic 
phase diagram (see figure caption for detailed parameters). We find that a highly-connected 
(large Pc) tense (short Pg) network consisting of stiff bonds (large /37) is located deeply 
inside the ferromagnetic phase region, implying that the model network would experience a 
true phase transition to a state with infinite correlation length and divergent relaxation time 
when the mean field configurational entropy vanishes. The associated ideal glass transition 
may underlie the observed glassy dynamics of living cells^. 

Moreover, the analogous magnets are located close to the strong glass forming liquid 
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FIG. 6. Locating analogous magnets on the RG phase diagram for a random field Ising model 
at zero average field (adapted from Fig. 3 in Ref. [3]). The colored symbols correspond to model 
networks with Lg = 1.8, Pc — 0.5 at /37 = 5 (green) and 20 (purple), and to those with Lg = 
1.2, /37 = 2.5 at Pc = 0.5 (blue), 0.8 (yellow) and 1 (red). The analogous magnets are located close 
to the strong glass forming liquid Ge02 (square mark). 

Ge02. Recent experiments on osmotically compressed cells by Weitz group^ suggest that 
cells under compressive stress indeed behave as strong colloidal glass formers, based on an 
observed Arrhenius-type exponential growth of viscosity for cells with increasing volume 
fraction, as well as their lower fragility compared to hard spheres. 



2. Variation of surface tension and interface width with reducing packing fraction: existence 
of a spinodal point 

(} The surface free energy density a 

In the same way that Cammarota et al^ computed the surface energy for hybrid inher- 
ent structures at zero temperature, we can define the surface free energy (Fg) as the free 
energy cost of forming a stable interface between mobile and localized regions, with respect 
to the total bulk energies of the two regions if they were to exist separately having homoge- 
neous low and high a values respectively. Mathematically, we define Fg = Ftot — {F^ + F^), 
where Fl^^^ = E.eL(^) [(3/2) ln(4(^V^) + E,6L(r) /?K//(|rt - rf |; ))] represents the 
bulk energy of the left (right) region with homogeneous {«J} ({ttt}) solutions (see the 
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initial configurations in Fig. [Tb), and Ftot denotes the total free energy of the stable inho- 
mogeneous a configuration in which there is a smooth interface interpolating between the 
pinned boundary layers (see the final configurations in Fig. Ub)- At finite temperatures, the 
stable a configurations represent a compromise between the cost of localizing a particle and 
the free energy gain realized by particles being able to avoid each other once localized. 




FIG. 7. The variation of interface energetics and structure with reducing packing fraction 77. (a) 
Initial and final a profiles along x axis at various packing fractions. Homogeneous high a and 
low a solutions get closer as 77 decreases, (b) Corresponding color-coded a configurations. Same 
color scheme as in Fig. [3l (c) Surface tension versus iteration round number. Final surface tension 
almost vanishes as 77 approaches 77^ 2^ 0.85. (d) Surface tension versus packing fraction. The model 
parameters are Lg = 1.2, = 2.5, Pc = 0.5. 

To investigate the possible existence of a spinodal point, we studied what happens as 
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we reduce the packing fraction by shrinking the particles while maintaining their number 
density. Another extreme case of great interest is the jamming limit where particles run 
into each other and get firmly stuck. In this case the reconfiguration barrier becomes too 
high for remodeling to be relevant on laboratory time scales. Thus in the jamming regime 
mean-field theories should be sufficient for many practical purposes. 

We show in Fig. [7] the initial and final a profiles at various packing fractions (panel a) and 
the corresponding color-coded a configurations (panel b). Panel d displays the surface free 
energy density (or surface tension) a defined as a = Fg/N^ obtained by SCP calculations (i.e. 
final values in panel c), versus the packing fraction rj. It is clearly seen that surface tension 
decreases with reducing packing fraction and ultimately vanishes at 77 c::^ 0.85. Concomitantly 
we find that, for r] < 0.85, high-a solutions corresponding to localized amorphous packing 
become mechanically unstable to even small thermal fiuctuations. This limit of stability for 
an amorphous system is comparable to the Lindemann criterion for melting of crystalline 
solids. Therefore, surface tension and mechanical stability of amorphous packing seem to 
vanish at a similar 77, the dynamic transition packing fraction ry^, which can therefore be 
identified as a spinodal point. Consistent with this interpretation we find, homogeneous 
high and low a solutions tend to merge as r] decreases toward tja- 

^ Quantifying the interface width ^ 

In Fig. [8](a) we plot log{a/aH) versus x along which spatial variation of a occurs. We 
denote by an (c^l) the high (low) a value at the pinned boundary to the left (right). Much 
above tja (see 77 = 1 case), the interface is quite sharp; in a single layer a changes from near 
an to near aL. Close to ija (see r] = 0.87 case), the transition becomes smoother with a 
slowly varying over several atomic layers. In both cases, there arises a surface tension a 
reflecting the deviation of a in the interfacial layers from the a for either of the bulk free 
energy minima. 

To quantify the interface width, we locate the interface at the position Xin (red vertical 
lines in Fig. EH) where the most significant spatial variation of a takes place. We then make 
a linear fit (black dotted lines) at Xin and define the interface width ^ as the inverse slope 
of the curve log(a/a/f) versus x, i.e. = "^^^(1/^^)^^ • As shown in Fig. [8](b), the 

X — Xifi 

interface gradually spreads out as the packing fraction decreases, corresponding to a lower 
free energy cost for a smooth spatially varying order parameter. At the spinodal point, a 
vanishing surface tension indicates the formation of a sufl&ciently smooth interface between 
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FIG. 8. The variation of interface width (^) with reducing packing fraction 77. (a) Quantification 
of the interface width at 77 = 1 (upper) and 77 = 0.87 (lower). Xin indicates the location of interface 
where most significant variation in a occurs. The black dotted line represents a linear fit to the 
\og(a/aH) versus x curve at Xin- (b) The interface width ^ in units of inter-particle spacing tq 
versus packing fraction. The interface broadens as 77 decreases toward 77^. Model parameters are 
Le = 1.2,/37 = 2.5,Pc = 0.5. 

bulk stable phases no longer costs any energy. 

Fig. M demonstrates how the glassy parameters vary with the packing fraction. As rj is 
reduced, configurational entropy density Sc increases (panel a) whereas surface tension J for 
single-particle droplets moderately decreases (panel b). But J does not vanish when a does. 
Instead, as 77 approaches the spinodal point (77^ = 0.85), J remains rather close to the Xia- 
Wolynes estimate. Moreover, the explicit mapping gives an entropy at crossover Sc{r]A) — 
2.8kB' This is approximately double the estimate given by the string transition theory^; the 
latter predicts barrierless reconfiguration events occur at a critical configurational entropy 
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^string ^ 1.13/^^. It is also quite a bit larger than the estimate based on percolation clusters 
^perc ^ l.28fc^. These results are consistent with the observation that activated dynamics 
and mode coupling effects coexist at the empirical crossover temperaturei^*^*^. 




FIG. 9. The variation of configurational entropy Sc (a) and surface tension J (b) with reducing 
packing fraction 77. The model parameters are Lg = 1.2, /37 = 2.5, Pc — 0.5. 



B. Active networks: dependence of glassy characteristics on motor properties 

While passive glassy systems remain trapped in micro-configurations for incredibly long 
times because only thermal energy is available to allow the system to evolve, ATP-dependent 
rearrangements of active cytoskeletal networks modify the micro-configurations themselves 
and provide an alternate means of exploring new configurations. ATP-hydrolysis-powered 
motor action can resolve steric constraints and drive structural rearrangements. 

1. Effect of motor action 

We have shown earlier^^ that a model cytoskeleton driven by small-step motor kicks can 
be described by an effective temperature Tgff that depends on motor activity A and motor 
susceptibility s through the relation T^f^/T = [1 + (5 — 1/2) A]""^. Using this explicit map- 
ping, we show in Fig. [10] that for a given A, as s decreases and thus Tgff rises, configurational 
entropy Sc increases (panels a and c) whereas the surface tension J decreases (panels b and 
d), indicating that the higher Teff induced by increasingly more load-resisting motors at a 
given activity promotes reconfiguration and reduces the mismatch free energy penalty. On 
the other hand, as seen in panel (b), for any given value of 5, higher motor activity leads 

16 



to a larger surface tension J . This is because intense motor action enhances the effective 
attraction in the buckhng regime, regardless of motor susceptibility. More interestingly, as 
shown in panel (a), while Sc decreases with increasing A for susceptible motors (5 > 0), 
load-resisting motors {s = —0.25) might resolve constraints and promote structural remod- 
eling, thus expanding the accessible configurational state space as motor activity increases, 
i.e. there is a larger Sc as A increases. 




Teff/T Xeff/T 



FIG. 10. Glassy characteristics for active networks. The dependence of configurational entropy 
density Sc (panels a and c) and surface tension J (panels b and d) on motor activity A, susceptibility 
s and effective temperature Teff is shown. We also indicate by black crosses the corresponding 
quantities for passive networks with the same underlying architecture (at A = and thus Teff = T) 
for comparison. Model parameters are Lg = 1.2, = 2.5, Pc = 0.5, 77 = 1. 

As a comparison, the black crosses in the figures show the values of Sc and J in the absence 
of motor processes, i.e. at A = and thus Teff = T. These correspond with what should 
be found after complete ATP depletion in experiments. We might say after "rigor mortis 
has set in" . The convergence of different curves to the passive values as A ^ is apparent 
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(see panels a and b). While panels (a) and (b) disentangle the effect of A and s on glassy 
dynamics, panels (c) and (d) demonstrate the use of Tgff in describing the non-equilibrium 
behavior driven by small-step motors: Raising Teff allows the system to explore a larger 
portion of the configuration space (larger Sc) and reduces the coupling strength (lower J). 



2. Estimate of the activation harrier 

RFOT theory allows us to estimate the magnitude of activation barriers. Using the 
expression from RFOT theory that includes the interface disorder wetting effect^ we can 
write 



log(r/ro) = /3effAF^ = 3^ 



(7) 



Fig. [TT] shows the dependence of the resulting activation barrier height /3effAF* on motor ac- 
tivity A (panel a) and configurational entropy Sc (panel b) for a series of motor susceptibility 
s. For susceptible motors with 5 > 0, more intense motor action raises the reconfiguration 
barrier. Load-resisting motors {s = —0.25) at sufficiently high activity, however, actually 
facilitate structural rearrangements thus lowering the barrier (see lowest curve in panel a). 
An inverse dependence of barrier height on the configurational entropy remains valid (see 
panel b). 



(a) 



o 




(b) 
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FIG. 11. The RFOT estimate of the magnitude of activation barrier. We plot barrier height 
log(T/To) versus (a) motor activity A and (b) configurational entropy Sc, for a series of motor 
susceptibility s. "Control" gives the passive value at A = for comparison. Model parameters are 
Le = 1.2, /37 = 2.5, Pe = 0.5, 77 = 1. 
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To translate the barrier height into relaxation times r in the laboratory, we need to 
estimate the bare relaxation time tq. tq = d\/ D^^^ is the time to diffuse a Lindemann length 
rf^, which is related to the localization strength an of confined motion through d\ = l/an- 
The eflFective diffusion constant given by D^^ = /c^T(l + A/3)/67r77or depends inversely on 
the viscosity r]o of the medium in which the network is immersed as well as the effective 
radius r of the crosslinks/nodes of the network. The drastic difference in relaxation time 
between molecular glass formers and colloidal glass formers lies in the big difference in their 

The Weitz group^ has recently shown that cells under compressive stress behave as strong 
colloidal glass formers; ATP depletion raises the reconfiguration barrier, whereas F-actin 
disruption accelerates remodeling. Our theory predicts a trend consistent with these obser- 
vations for the dependence of glassy dynamics on the nonequilibrium processes and network 
connectivity. Fig. [12] demonstrates how the activation barrier height varies with packing frac- 
tion at various network connectivities and motor susceptibilities. The passive case (black, 
circle) without motor action corresponds to ATP depletion in experiment. Compared to 
the active network at the same connectivity {Pc = 0.5) yet driven by load-resisting motors 
{s = —0.25; blue square), the lack of motor processes in the passive case indeed raises the 
activation barrier. Such motor-facilitated remodeling is also consistent with active fiuidiza- 
tion observed in polymer networks where myosin-II motors enhance longitudinal filament 
motion^^. On the other hand, a reduction in connectivity {Pc = 0.4, s = —0.25; red square) 
mimicking F-actin disruption leads to a faster relaxation, since fewer bond constraints yield 
weaker interactions and larger configurational entropies and hence a lower reconfiguration 
barrier. 

In contrast to the load-resisting motors, at a given connectivity, susceptible motors = 1, 
lines with triangular marks) raise the activation barrier and suppress remodeling. This 
confiicting behavior suggests a way to determine the sign of the effective motor susceptibility 
in living cells. The experimental observation that active processes promote reconfiguration 
for a broad range of cell types^ thus suggests molecular motors in the cytoskeleton tend to 
work against mechanical load {s < 0) rather than move energetically downhill {s > 0). As 
we have shown earlier, load-resisting motor action is also essential for macroscopic active 
contractility in actomyosin networks^^. Finally, in all the cases, reconfiguration slows down 
(i.e. r increases) as the volume fraction rj increases, as seen experimentally^. 

19 



6 



p^=0.5, s=l 



5 - 



p^=0.5, s=-0.25 
p^=0.4, s=-0.25 




1 - 



0.85 



0.9 



0.95 



packing fraction r| 



FIG. 12. Dependence of reconfiguration barrier on motor susceptibility and network connectivity. 
We plot barrier height log(T/To) versus packing fraction 77 for susceptible motors (5 = 1, triangle) 
and load-resisting motors (s = —0.25, square) at different network connectivity Pc. Passive curve 
(with black circles) corresponds to ATP-depletion in experiments. Model parameters are Lg = 
1.2, /37 = 2.5, A = 1. 

VI. SUMMARY 

In this work, we have developed a theoretical framework to study the glassy dynamical 
behavior of passive and active network materials. To illustrate the concepts, we have built a 
general microscopic network model that incorporates both the nonlinear elasticity of individ- 
ual filaments and the steric constraints due to crowding. By treating the network materials 
as structural glass forming liquids, we draw an explicit analogy to disordered short-range 
Ising magnets. This magnetic analogy, when combined with renormalization group phase di- 
agrams for a random-field Ising model, allows us to predict whether a thermodynamic phase 
transition underlies the observed glassy dynamics in network materials. The calculations 
also locate the model networks close to strong glass formers, as found experimentally. 

Moreover, this explicit mapping provides a microscopic route to compute the mismatch 
surface tension of RFOT theory and the configurational entropy and thus lends insight on 
how the network architecture and motor properties influence the activation dynamics asso- 
ciated with structural rearrangements. Consistent with experimental observations in living 
cells, we find that intense action of load-resisting motors allows the systems to explore a 
larger portion of the configurational space and thus facilitates more efficient remodeling. A 
heterogeneous self-consistent phonon approximation enables us to study spatially varying 
order parameters with a smooth gradient across the system. This method proves useful in 
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identifying a spinodal point where formation of a sufficiently smooth interface between local- 
ized and freely moving particles no longer costs any energy and activated dynamics crosses 
over to collisional transport. Finally, near-universality of the interaction per molecular unit 
based on universality of Lindemann ratio of all glass formers made of spherical particles still 
proves to be a good approximation in deeply glassy regime (i.e. at sufficiently high packing 
fraction) for both passive and active network materials. 
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